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We show that lateral fluidity in membranes containing quenched protein obstacles belongs to the 
universality class of the two-dimensional random-field Ising model. The main feature of this class is 
the absence of a phase transition: there is no critical point, and macroscopic domain formation does 
not occur. Instead, there is only one phase. This phase is highly heterogeneous, with a structure 
consisting of micro-domains. The presence of quenched protein obstacles thus provides a mechanism 
to stabilize lipid rafts in equilibrium. Crucial for two-dimensional random-field Ising universality is 
that the obstacles are randomly distributed, and have a preferred affinity to one of the lipid species. 
When these conditions are not met, standard Ising or diluted Ising universality apply. In these 
cases, a critical point does exist, marking the onset toward macroscopic demixing. 

PACS numbers: 87.16.dt, 64.60.-i, 61.20.Ja 



I. INTRODUCTION 

Membranes are two-dimensional (2D) fluid environ- 
ments pQ. There is growing consensus that membrane 
lateral structure is heterogeneous, and characterized by 
domains of different size and composition [2-6 . Domain 
formation in biological membranes is important because 
it links to key processes in cells, such as signaling, endocy- 
tosis, and adhesion [3 [8], while in model membranes do- 
main formation is relevant for applications, ranging from 
photolithographic patterning, spatial addressing, micro- 
contact printing, and microfluidic patterning [9] [10] . To 
identify the factors that control domain formation, and to 
understand the underlying physical mechanisms, is there- 
fore of practical importance. 

In thermal equilibrium, a heterogeneous structure is 
difficult to comprehend, due to the large cost in line ten- 
sion [2 j. One instance where a heterogeneous structure 
does arise is near a critical point, as was demonstrated 
in free-floating giant unilamellar vesicles (GUVs) [TTHT3] . 
Near a critical point, the line tension vanishes; thermal 
motions then induce composition fluctuations over a wide 
range of length scales. An important concept in the the- 
ory of critical phenomena [T3j [14] is universality: systems 
that belong to the same universality class undergo sim- 
ilar phase transitions, and yield the same set of critical 
exponents. For GUVs, the universality class was shown 
to be that of the 2D Ising model [TTj . 

Since the universality class does not depend on the mi- 
croscopic details of a system, it is tempting to speculate 
that 2D Ising universality, as observed in GUVs, is the 
generic class for membrane fluidity. The purpose of this 
paper is to show that, in less-idealized membranes (com- 
pared to GUVs), a different universality class comes into 
play, the reason being the presence of quenched obsta- 
cles. In living cells, proteins can bind to the underly- 
ing cytoskeleton [15 . This leads to a 2D fluid consist- 
ing of mobile particles (e.g. lipids) diffusing in a back- 
ground of quenched (immobilized) protein obstacles. A 



similar situation arises in supported membranes, where 
surface friction may lead to particle immobilization [16] . 
Hence, to understand domain formation in membranes, it 
is important to keep in mind that an immobilized com- 
ponent may be present (in physical terms, such a sys- 
tem is called a quenched-annealed mixture [17]). Pre- 
cisely this point was recognized by two recent simulation 
studies [18] [19] where domain formation in membranes 
with quenched protein obstacles was investigated. It was 
found that quenched obstacles lower the critical temper- 
ature [T8l[T9]. 

The aim of this paper is to relate the findings of 
Refs. 18] and [19] to universality classes. Our main mes- 
sage is that, in the presence of quenched obstacles, mem- 
brane fluidity belongs to the universality class of the 
2D random-field Ising model (2D-RFIM). The 2D-RFIM 
is crucially different from the 2D Ising model because 
it does not feature a phase transition [20] [21] . Conse- 
quently, critical behavior and macroscopic phase sepa- 
ration do not occur, even in the limit of low obstacle 
concentration. This means that, irrespective of temper- 
ature and lipid composition, a membrane with quenched 
obstacles is always in the same thermodynamic phase. 
The structure of this phase is found to be heterogeneous, 
consisting of micro-domains. Hence, based on univer- 
sality alone, one elegantly accounts for a heterogeneous 
equilibrium domain structure over a wide range of com- 
positions and temperatures. In contrast, critical fluctu- 
ations persist only in a small region around the critical 
point. 

The universality class of the 2D-RFIM applies when 
the protein obstacles display a preferred affinity to one 
of the lipid phases, and are randomly distributed. We be- 
lieve this to be the typical situation in most membranes. 
When these conditions are not met, standard 2D Ising 
(as observed in GUVs) or diluted 2D Ising universality 
[22] will arise. However, deviations between the latter 
two classes are relatively small, in the sense that both 
feature a critical point, with similar critical exponents 
[231 [24]. 
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II. MODEL AND METHODS 

To illustrate these points, computer simulations of a 
2D off-lattice mixture with quenched obstacles are per- 
formed. Since we focus on universality - which does not 
depend on microscopic details - it suffices to use a sim- 
ple model. The protein obstacles are unit diameter disks, 
and are placed at the start of each simulation run on a 
L x L square with periodic boundaries, after which they 
remain quenched. The obstacles are mostly distributed 
randomly, although non-random choices will also be in- 
vestigated. In line with recent simulations of membrane 
criticality [TTHT31 [25] a two-state model is used to de- 
scribe the lipids, which diffuse through the environment 
of quenched obstacles. The model contains saturated (S) 
and unsaturated (U) lipids, both of which are disks hav- 
ing the same diameter as the obstacles. The sole inter- 
action between lipids is a hard-core repulsion between 
unlike species, which is a minimum condition to induce 
phase separation [26^ In addition, the lipids interact 
with the quenched obstacles; these interactions will be 
specified later. 

In ternary mixtures of saturated and unsaturated lipids 
with cholesterol, critical behavior occurs over a wide 
range of compositions and temperatures [27 . To in- 
duce critical behavior, one either varies the temperature 
at fixed composition, or one varies the composition at 
fixed temperature; both routes are experimentally acces- 
sible [28 . Since the interaction in our model is hard- 
core, we vary the lipid composition. This is done in 
grand-canonical Monte Carlo (MC) simulations, i.e. the 
respective fugacities zs and zjj of saturated and unsat- 
urated lipids are fixed, but the lipid number densities 
px = Nx/L 2 fluctuate, with Nx the number of lipids 
of species X G (U, S) in the system (see Ref. [29l for full 
details on the simulation method). The phase behav- 
ior is analyzed using the order parameter distribution 
(OPD) Pl,i(pu\zui z s), defined as the probability to ob- 
serve a state with unsaturated lipid density pu- The 
OPD depends on the imposed fugacities, the system size 
L, as well as on the configuration i of quenched obsta- 
cles. Since universality applies to the thermodynamic 
limit L — >• oo, meaningful results require simulation data 
over a range of L, and the use of finite-size scaling (FSS) 
[30] to perform the extrapolation L — >• oo. In addition, 
in the presence of quenched obstacles, the shape of the 
OPD may fluctuate profoundly between obstacle config- 
urations. It is therefore crucial to average simulation 
results over many i — 1, . . . , K different obstacle config- 
urations. 



III. RESULTS 

A. pure lipid membrane: 2D Ising universality 

We first consider a pure lipid membrane, i.e. without 
quenched obstacles (this qualitatively resembles a GUV) . 
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FIG. 1. Finite-size effects in a pure lipid membrane at fu- 
gacity zs — 2 where the demixing transition is first-order, 
(a) The OPD for L = 50, and with zu tuned according to 
Eq. ([!]). The distribution is distinctly bimodal. The vertical 
arrow marks the free energy barrier AFl of interface forma- 
tion, (b) The variation of AFl with L. The linear increase 
confirms Eq.|2| of a genuine first-order transition; from the 
slope of the line the line tension can be extracted. The lower 
frames show snapshots of the membrane corresponding to the 
left peak of the OPD (c), the coexistence region between the 
peaks (d), and the right peak (e), where dark (light) particles 
represent saturated (unsaturated) lipids. 



We thus simulate saturated and unsaturated lipids only, 
both species being mobile. The control parameters are 
the particle fugacities, zs and zjj, and by appropri- 
ately tuning these parameters, phase transitions (pro- 
vided they exist) can be induced. In some sense, the pure 
lipid membrane may be conceived as an off-lattice version 
of the 2D Ising model, with the role of temperature re- 
placed by 1/zst an d that of the external field by \rizu- 
When zs exceeds a critical value zs > z CVl the tran- 
sition is first-order and macroscopic phase separation is 
observed. When zs = z CTl the transition becomes contin- 



3 



1.5 



L=20 
L=25 
L=30 
L=35 
L=40 
L=45 
L=50 
2D Ising value 




FIG. 2. Demonstration of the cumulant intersection method 
[311 [32] to locate the critical point in a membrane without 
quenched obstacles. Plotted is the cumulant Ui versus the 
saturated lipid fugacity zs for several system sizes L, and zjj 
tuned according to Eq.|l]). At the critical point, the curves 
for different L intersect. The horizontal line marks U\ ~ 1.07 
of the 2D Ising universality class. 



uous, and critical behavior is observed. When zs < z CT1 
there is only one phase. At a phase transition, the scal- 
ing of the OPD with L assumes a characteristic form, 
from which the transition type can be determined. The 
observation of a transition also requires that zjj is chosen 
suitably. In what follows, for a given value of z$ : we tune 
zjj such that the derivative of the average saturated lipid 
density with respect to zjj is maximized 



d(pu) s 
zjj : tt, > max, 

am Zjj 



(1) 



as is commonly done in fluid simulations [33 . In the 
above, (•) is a thermal average, i.e. an integral over the 
OPD: (pjj) = J °° pu P L (p u \z u ,z s )dp u . 

At a first-order transition the OPD is bimodal; see 
Fig.[TJa), where the natural logarithm of the distribution 
is shown (there is obviously no dependence on the obsta- 
cle configuration i in this case). Provided the bimodal 
shape persists in the thermodynamic limit L — >• oo, the 
peaks reflect stable phases. The peak on the right cor- 
responds to a homogeneous phase rich in unsaturated 
lipids, and lean in saturated lipids (Fig. [TJc)). The left 
peak corresponds to a homogeneous phase of reversed 
composition (Fig. [TJe)), i.e. rich in saturated lipids, and 
thus resembles a lipid raft [5]|34[|35]. However, rafts are 
defined to be micro-domains, while the peaks in the OPD 
reflect macroscopic phases. 

To verify that the bimodal shape survives in the ther- 
modynamic limit, we consider the variation of the peak 
height AFl with L (vertical arrow in Fig. [TJa)). We em- 
phasize that AFl is obtained from the natural logarithm 
of the OPD: it is defined as the average peak height, mea- 
sured from the minimum "in-between" the peaks. When 
the simulation traverses the region between the peaks, 



phase coexistence is observed. Both phases then appear 
simultaneously. Exactly between the peaks, each phase 
occupies half the system: the phases then arrange in two 
slabs since this yields the shortest interface (Fig. [TJd)). 
Note that, due to periodic boundaries, two interfaces are 
present. If L is large enough, the interfaces do not inter- 
act with each other: the relative amount of the phases 
can then be varied over some range without any cost in 
free energy, which is the origin of the characteristic flat 
region between the peaks in the distribution of Fig. [TJa) 
[36] . Provided a flat region between the peaks is present, 
AFl corresponds to the free energy cost of interface for- 
mation [37] . Since, in the slab arrangement, the total 
interface length / s i a b = 2L, it follows that 



AF L = T^siab 27L, 



(2) 



with 7 the line tension [37] [38] . At a first-order transi- 
tion in 2D, we thus expect a linear increase of AFl with 
system size. This is confirmed in Fig.[TJb), and from the 
slope of the line 7 can be determined. Note that, for 
small L, corrections to Eq.Q become important [39] , as 
indicated by the systematic deviation of the simulation 
data away from a straight line for L < 25 or so. 

Precisely at zs = z CTl the first-order transition termi- 
nates in a critical point. The hallmark of criticality is 
scale invariance. This property can be exploited in FSS 
to extract z cr from simulation data [3TJ[32]. To this end, 
one measures the Binder cumulant U\ = (ra 2 )/(|m|) 2 , 
m = pu — (pu), which in the thermodynamic limit as- 
sumes three distinct values [3D]. For zs > z cvi i.e. where 
the transition is first-order and characterized by two- 
phase coexistence, the OPD is bimodal (Fig. [TJa)). The 
OPD may then be approximated by a superposition of 
two non-overlapping Gaussian peaks [40], which can be 
shown to yield Ui = 1. For zs < z CTl there is only 
one phase, with the OPD, consequently, featuring just 
a single peak. The peak is again Gaussian, for which 
U\ = 7r/2. Precisely at the critical point, the OPD is 
bimodal, but the peaks overlap. As a result, the Binder 
cumulant at criticality assumes an "in-between" value 
1 < Ui < 7r/2; scale invariance implies that [/* does 
not depend on L. Moreover, U* is characteristic of the 
universality class, and for the 2D Ising model Ui ~ 1.07. 

In the thermodynamic limit, the Binder cumulant thus 
equals: 



lim Ui 

L— >oo 



1 Z S > Z c 

Ui z s z c 
7r/2 z s < z c 



(two-phase region), 
(critical point), (3) 
(one-phase region). 



This behavior is well-suited to obtain the critical fugacity 
z cr from simulation data. To this end, one plots U\ versus 
zs for different system sizes (Fig.|2|. At the critical point, 
the curves for different L intersect. From the intersec- 
tion point we conclude that z CY ^1.72 for the membrane 
model without quenched obstacles [29]. Note also that 
Ui at the intersection is close to the 2D Ising value (hor- 
izontal line). Hence, from the Binder cumulant alone, we 



(a) (b) 
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FIG. 3. Some example OPDs for a membrane with quenched 
obstacles, with the obstacles randomly distributed, and with 
a preferred affinity to saturated lipids. The key observation is 
that the distributions profoundly fluctuate between obstacle 
configurations, which illustrates the need for an extensive av- 
erage over many obstacle configurations; data are shown for 
L = 20, zs — 1.92, and zjj tuned according to Eq.|l]). 
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FIG. 4. Cumulant analysis of a membrane with quenched 
obstacles, with the obstacles randomly distributed, and with 
a preferred affinity to saturated lipids. Plotted is U[ versus 
zs for different L. An intersection point is not observed: this 
indicates the absence of a phase transition, which is consistent 
with the universality class of the 2D-RFIM. 



prove the existence of a critical point, obtain an estimate 
of the critical fugacity, and confirm 2D Ising universality. 
The critical fugacity, as well as the total lipid density at 
criticality p cr ^1.56 (which was extracted from the criti- 
cal OPD), are in good agreement with other studies [SI] . 
Of course, since the pure membrane model is symmetric 
under the inversion of U S lipids, it trivially follows 
that ps = pu at the critical point. 

B. random obstacles with preferred affinity: 
2D random-field Ising universality 

We now come to the main result of this paper, where 
we consider a membrane with quenched protein obsta- 
cles. The obstacles have a hard-core interaction with the 
unsaturated lipids, while saturated lipids may overlap 
freely with them. We thus impose a preferred affinity of 
the obstacles to lipid rafts, as, for example, GPI anchored 
proteins are known to do [34] . At the start of each simula- 
tion, we place Q obstacles randomly in the system, after 
which they remain quenched; Q is drawn from a Pois- 
son distribution P(Q) = (\L 2 )®e- XL2 /Q\ with A = 0.03 
the average obstacle density. Next, saturated and unsat- 
urated lipids are introduced, and the OPD is measured 
for the given configuration of quenched proteins. This 
procedure is repeated for many different obstacle config- 
urations. 

Since the proteins have a preferred affinity to one of the 
lipid species, and since they are randomly distributed, 
we expect 2D-RFIM universality. In particular: there 
should no longer be a phase transition. In Fig. [3] OPDs 
for a number of obstacle configurations are shown. The 
distributions were obtained for zs = 1.92, which signif- 



icantly exceeds z CY of the membrane without quenched 
obstacles. We observe large shape variations. In (a) we 
see a bimodal distribution, which might be taken as ev- 
idence of a first-order transition. However, in (b) only a 
single peak is revealed, which rather reflects no transition 
at all. Finally, (c) and (d) show distributions with three 
peaks, suggesting a triple-point. We emphasize that all 
distributions in Fig. [3] were obtained at the same "inverse 
temperature" zs and system size: only the obstacle con- 
figurations are different. Since the distributions show ex- 
treme shape variations between obstacle configurations, 
it is clear that meaningful results require an average over 
many obstacle configurations. One might object that the 
shape variations in Fig. [3] merely reflect a finite-size ar- 
tifact, and that an average over obstacle configurations 
is not needed in larger systems. This, however, is a dan- 
gerous assumption because random-field systems are no- 
toriously non- self- averaging [42] [43] . 

We thus use K = 2000 obstacle configurations in what 
follows, and perform a FSS analysis to obtain insight 
in the phase behavior in the thermodynamic limit. To 
this end, we consider the quenched- averaged cumulant 
U[ = [(^ 2 )]/[(|^|) 2 ], with m and (•) defined as be- 
fore, and where [•] is an average over obstacle config- 
urations. For each obstacle configuration z, the OPD 
is tuned according to Eq.Q [44 , which is then used 
to compute the thermal averages (•)$; the latter are 
subsequently averaged over the obstacle configurations 

[(')] = (lfflEf=i(i In Fig. Q we plot U[ versus z s 
for different L. In contrast to Fig. [2] we cannot identify 
an intersection point. Instead, for a fixed value of zs, the 
trend is that, by increasing L, the cumulant approaches 
the value of the one-phase region: lim^^oo U[ — » tt/2. 
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FIG. 5. Computer generated snapshots of domain formation 
in a membrane with quenched obstacles; L = 60, p = 1.6, and 
A = 0.04 were used. The dark (light) regions correspond to 
saturated (unsaturated) lipids; dots mark the protein obsta- 
cles. 



Hence, irrespective of zs, in the thermodynamic limit 
L — >> oo there is only one phase. In agreement with 2D- 
RFIM universality, a phase transition (above which co- 
existence between two macroscopic phases would occur) 
no longer takes place. 

The reason macroscopic phase separation in the 2D- 
RFIM is prevented has a clear physical origin: in the 
presence of quenched obstacles the cost of line tension 
becomes negligible. To see this, consider a raft domain 
of radius R. The cost of line tension scales ex R d ~ x , 
with d the spatial dimension. However, the raft will also 
encompass a number of obstacles; the typical number 
of obstacles in the domain is oc R d but with Poissonian 
fluctuations oc R d l 2 . Hence, it is favorable for the raft 
to seek out those regions in the membrane where there 
is an excess of obstacles. The cost of line tension is then 
compensated by the quenched obstacle excess, since d = 
2 for membrane fluidity. Note that this reasoning is just 
the analogue of the Imry-Ma argument for random-field 
magnets [45] . 

Having shown that randomly distributed quenched 
proteins (with preferred affinity to one of the lipid 
species) prevent macroscopic phase separation, the struc- 
ture of the remaining single-phase will now be analyzed. 
To this end, we perform canonical (fixed density) MC 
simulations using ps = pu f° r different total densities 
P : = Ps + Pu- In Fig. [5] we show a number of equilibrated 
snapshots at p = 1.6, L = 60, and using the same obsta- 
cle configuration each time; the obstacle concentration 
A = 0.04. Note that p = 1.6 exceeds the critical density 
of the pure membrane and so, if we were to remove the 
obstacles, the snapshots would reveal macroscopic phase 
separation, and resemble Fig. [TJd). However, from the 
cumulant analysis of Fig. [3J it is clear that this does not 
happen in the presence of quenched proteins: there is no 



longer a critical point, and hence no longer a two-phase 
coexistence region. Instead, we see a structure consist- 
ing of micro-domains (by micro we mean that the typical 
domain size exceeds the typical distance between the ob- 
stacles but remains finite). We also see that the domain 
structure in each of the snapshots is different: this shows 
that thermal fluctuations are still present, and that the 
domains have a finite lifetime (which is important be- 
cause rafts are believed to be short-lived). At first sight, 
the domains in the snapshots of Fig. [5] look deceivingly 
similar to critical fluctuations; see for example Fig. 1 of 
Ref. |TT| where critical fluctuations in GUVs are shown. 
However, the domains in Fig. [5] are crucially different in 
two respects. 

The first difference is that critical fluctuations are spa- 
tially indifferent: they can form at any location in the 
membrane with equal probability. In contrast, the do- 
mains that form in the presence of quenched obstacles are 
spatially selective: following the Imry-Ma argument, raft 
domains prefer to form at those locations in the mem- 
brane that feature an excess of obstacles. Of course, 
this spatial selection does not appear directly in single 
snapshots, but it becomes strikingly visible when we con- 
sider many snapshots and average over them. To demon- 
strate this explicitly, we use the same obstacle configura- 
tion of Fig. [5j and create a large number of equilibrated 
snapshots. For this set of snapshots, we construct an 
"ensemble- averaged" snapshot, simply by overlaying the 
individual snapshots. To this end, a grid consisting of 
unit square cells is placed over each individual snapshot. 
For each grid cell /c, we sum over all snapshots, and count 
how often the cell contained a saturated lipid, or an un- 
saturated lipid; the counts are denoted Cs(k) and Cu(k), 
respectively. From these counts, we compute the canon- 
ical (or time) averaged lipid preference for each cell 



A k := 



Cu(k)-C s {k) 

Cu(k) + c s (ky 



(4) 



where a value Ak — — 1 (Ak — 1) indicates a preference 
of cell k to saturated (unsaturated) lipids, while Ak = 
indicates that a preferred affinity is absent. 

In case of critical fluctuations, which are spatially in- 
different, the structure revealed in individual snapshots 
gets completely "washed-out" in the averaged snapshot, 
since Ak = for all cells in that case. In contrast, for 
the membrane with quenched obstacles, a clear structure 
in the averaged snapshot remains visible, which marks 
the regions in the membrane where rafts are most likely 
to be found. In Fig. |6j we show a number of averaged 
snapshots thus obtained, each one using the same obsta- 
cle configuration, but at different lipid densities p. The 
snapshots are color-coded, and the color reflects the value 
of Ak in the grid cell; dots mark the quenched protein 
obstacles. For all lipid densities considered, which even 
includes one density slightly below the critical density of 
the pure model, it is clear that lipid domains are spatially 
selective. 

The second difference is that critical fluctuations van- 
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FIG. 6. Series of ensemble- averaged snapshots of a membrane with quenched protein obstacles, for three values of the lipid 
density p, using L = 60 and A = 0.04. The images are color-coded, and the color reflects the local lipid preference Ak defined 
in Eq.Q, from Ak = — 1 (brightest) to Ak = 1 (darkest); dots mark the protein obstacles. Note that the same obstacle 
configuration is used in all three cases; for a different analysis where also an average over obstacle configurations is performed, 
see Fig. [7] The images clearly show that the proteins create regions in the membrane where certain lipid species are preferred. 
In addition, the shape of these "preferred" regions is remarkably stable under variations in the lipid density. 



ish above the critical density: the heterogeneous domain 
structure is then replaced by a coexistence between two 
(homogeneous) macro-domains. In contrast, the micro- 
domain structure in the membrane with quenched ob- 
stacles is stable upon variations in the total lipid den- 
sity. The main effect of increasing p is a "freezing-out" of 
thermal fluctuations. At high density, the membrane be- 
comes more spatially selective, and it becomes less likely 
to observe raft domains in regions where there is no ob- 
stacle excess. This is manifested in the averaged snap- 
shots of Fig. [6] by a "sharpening" of the domain walls. 

This effect can be made more precise when we convert 
each averaged snapshot of Fig. [6] into a one-dimensional 
histogram of Ak values, and to subsequently average 
these histograms over a number of obstacle configura- 
tions. Of course, in this way we loose the spatial cor- 
relations, but it allows us to verify that the spatial se- 
lectivity shown in Fig. [6] is a generic feature, and not 
just an artifact of the particular obstacle configuration 
that was used. In Fig. [7| we show histograms of Ak 
values thus obtained, which were averaged over 20 dif- 
ferent protein configurations, and again for three lipid 
densities p. Note that the sharp peak at Ak = —1, in 
particular for the p = 1.4 histogram, is partially due to 
grid cells containing an obstacle. By increasing p, the 
extreme values Ak ~ ±1 become more likely, which con- 
firms that the preference for certain lipids at certain loca- 
tions in the membrane becomes more pronounced. The 
histograms also reveal another important point: even at 
density p = 1.4, i.e. below the critical density of the pure 
model, the membrane with quenched proteins is already 
spatially selective, since no peak around Ak = is visible. 

To summarize: in the presence of quenched proteins, 
randomly distributed, and with a preferred affinity to one 
of the lipid species, there is no sign of a phase transition 
leading to two-phase coexistence, nor of a critical point. 
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FIG. 7. Histograms of observed Ak values for different values 
of the lipid density p collected over 20 obstacle configurations. 



Instead of the lipids forming two macroscopic domains, 
micro-domains are observed. These micro-domains are 
dynamic, but they do not diffuse over the membrane ran- 
domly: they are most likely to be found in regions of 
the membrane featuring an excess of obstacles. In some 
sense, the quenched protein obstacles provide a scaffold- 
ing (channels) extending over the membrane for raft do- 
mains. 
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FIG. 9. Variation of the critical fugacity z cr with obstacle 
concentration A for a membrane with neutral obstacles that 
are randomly distributed (squares), and for a membrane with 
obstacles featuring a preferred affinity to saturated lipids, but 
placed on a regular grid (triangles). The dot at A = marks 
z C r of the pure membrane; lines serve to guide the eye. 



FIG. 8. (a) Cumulant analysis of a membrane with quenched 
protein obstacles that are randomly distributed, but without 
a preferred affinity to one of the lipid species. The curves of U[ 
versus zs for different L reveal an intersection point, consis- 
tent with the occurrence of a phase transition, (b) The same 
analysis as above, but for protein obstacles with a preferred 
affinity to saturated lipids placed on a regular grid. Curves 
of Ui versus zs for different L intersect, consistent with the 
occurrence of a phase transition. 



C. random obstacles without preferred affinity : 
diluted 2D Ising universality 

We now consider a membrane with protein obstacles 
that do not have a preferred affinity to one of the lipid 
species, i.e. are neutral, but remain randomly distributed. 
To this end, we use the rule that neither lipid species may 
overlap with the obstacles; the average obstacle concen- 
tration A = 0.03, again drawn from a Poisson distribu- 
tion. The Imry-Ma argument does not apply in this case, 
and we indeed find a radical departure from 2D-RFIM 
universality. Instead, the scenario of Fig. [I] is recovered, 
and a genuine phase transition is observed. In Fig. [8ja), 
we show the variation of the cumulant U[ with zs, for 
different system sizes. An intersection point is revealed, 
proving the existence of a critical point, and for the criti- 
cal fugacity we obtain z cr w 1.86. Since there is no longer 
a preferred affinity, the protein obstacles cannot compen- 
sate the cost of line tension. The membrane thus seeks 
to minimize the amount of interface, implying the forma- 
tion of macroscopic domains when zs > z cr . Hence, in 



the presence of neutral obstacles, the only instance where 
a heterogeneous domain structure in equilibrium can arise 
is near the critical point. The corresponding critical fluc- 
tuations will be qualitatively similar to those of the pure 
membrane, i.e. they are dynamic, and spatially indiffer- 
ent (since the membrane with neutral obstacles remains 
symmetric under inversion of U «->• S lipids, it trivially 
follows that the averaged lipid preference Ak =0). 

We emphasize that the case of neutral obstacles closely 
resembles the membrane simulations of Ref. [18j In that 
work, the 2D Ising model is studied, but with a fraction 
of randomly chosen lattice sites "turned-off" to repre- 
sent quenched protein obstacles. The key point is that 
this kind of dilution disorder does not break the up/down 
symmetry of the Ising model, just as the neutral obsta- 
cles in our model do not break the U <H> S lipid sym- 
metry. Consequently, both models belong to the same 
universality class, namely the one of the 2D diluted Ising 
model. For this universality class, there is no doubt that 
a phase transition exists [22J, provided the dilution re- 
mains below the limit where the lattice becomes disjoint. 
Indeed, an analysis of the Binder cumulant in Ref. [18] 
also reveals an intersection point, in agreement with our 
Fig. (8^a). A second hallmark of diluted Ising universal- 
ity is a pronounced decrease of the critical temperature 
with increasing obstacle concentration. This was strik- 
ingly confirmed in Ref. HH and our data reveal the same 
trend (Fig. [9]). Note that zs in our model plays the role 
of inverse temperature, so z cr increases with the obstacle 
concentration. Finally, we point out that the extreme 
fluctuations of the OPD between obstacle configurations 
that characterize the random-field case (Fig. [3| do not 
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FIG. 10. Structure analysis of a membrane containing 
quenched protein obstacles; the obstacles have a preferred 
affinity to saturated lipids, and are placed on a grid with 
a lattice constant of five particle diameters. The left frame 
shows a single snapshot taken at densities slightly below the 
critical values pu,cr ~ 0.82 and ps,cr ~ 0.90 (due to the pre- 
ferred affinity, it does not hold that pu,cr = ps,cr)- The right 
frame shows the corresponding ensemble- averaged snapshot. 
The main message is that the structure seen on the left is 
completely washed out. All that remains is a locally enhanced 
saturated lipid density around each obstacle. 



occur in the diluted Ising model. 



D. non-random obstacles with preferred affinity : 
2D Ising universality 

Finally, we consider a membrane containing quenched 
proteins with a preferred affinity to saturated lipids, but 
placed non-randomly, namely on the sites of a square 
lattice. The lattice constant is six particle diameters, 
corresponding to an obstacle density A « 0.03. In this 
case, the Imry-Ma argument does not apply either, be- 
cause local (Poissonian) fluctuations in the obstacle den- 
sity do not occur; note also that a disorder average [•] is 
not needed either. In Fig. |8jb), we show the cumulant 
Ui versus zs for different L. We observe an intersection 
point, which proves the existence of a critical point, and 
for the critical fugacity z cr « 1.83 is obtained. Hence, the 
scenario of Fig. [T] is recovered, with macroscopic phase 
separation taking place when zs > z cr . By placing the 
protein obstacles at regular locations, the universality 
class remains that of the standard 2D Ising model. The 
only effect the proteins induce is confinement of the fluid, 
leading to a decrease of the critical temperature. Indeed, 
if we repeat the analysis using different lattice constants 
for the obstacle grid, we systematically find that z cr in- 
creases with obstacle concentration (Fig. [9]). 

Hence, with the obstacles placed on a grid, an equilib- 
rium heterogeneous domain structure only survives at the 
critical point. A typical snapshot obtained near critical- 



ity is given in the left frame of Fig. 10 , where now lattice 
constant five was used. Even though the obstacles have 
a preferred affinity to saturated lipids, this is not suffi- 
cient to stabilize a structure consisting of micro-domains. 



This becomes clear when one constructs the correspond- 
ing ensemble- averaged snapshot (conform Fig. |6| which 
we show in the right frame of Fig. [10] In the aver- 
aged snapshot, the structure seen in the single snapshot 
is completely washed out. Of course, due to the pre- 
ferred affinity to saturated lipids, there is a layer of en- 
hanced saturated lipid density around each obstacle, but 
this layer is restricted to a single obstacle. These layers 
do not extend to other obstacles to form micro-domains 
(one might regard them as nano-domains instead). In 
contrast, the micro-domains in Fig. [6] extend over nu- 
merous protein obstacles. 



E. Conclusion 

Inspired by recent simulations of membranes with 
quenched protein obstacles [I8j [19], additional simula- 
tions were performed with the aim to relate domain for- 
mation in these systems to universality classes. The main 
message is that, while a membrane without quenched ob- 
stacles is in the universality class of the 2D Ising model 
[TT] , the introduction of obstacles will in most practi- 
cal situations induce a change to 2D random-field Ising 
universality. The consequences are drastic: the phase 
transition of the pure membrane without obstacles is 
destroyed, macroscopic phase separation no longer oc- 
curs, and only a single phase survives consisting of micro- 
domains. These micro-domains prefer to form at special 
regions in the membrane, and this preference becomes 
more pronounced with increasing lipid density. It would 
be of great interest to see whether transport of particles 
happens along the same regions; in that case, quenched 
proteins truly provide a structure of channels along which 
membrane constituents are laterally transported. 

We have also considered two somewhat "artificial" 
cases, whereby the quenched obstacles are neutral, or 
placed regularly on a grid. In these cases, a normal 
critical phase transition is observed, above which macro- 
scopic phase separation takes place. These predictions 
are presumably less relevant for biological cells, but they 
could be verified in experiments on model membranes 
(using, for example, an appropriately patterned surface 
[46]). We write presumably above because the spectrin 
cytoskeleton network in red blood cells does, in fact, form 
a regular lattice structure [47] . 

On a more fundamental level, the connection between 
a fluid with quenched obstacles (with preferred affinity) 
and random-fields dates back to de Gennes [48] , but has 
been notoriously difficult to verify experimentally (see 
Ref . [49] for a review) . An alternative for the hypothe- 
sis of de Gennes, which in some cases better captures 
experiments, is the "wetting hypothesis" of Ref. [50] In 
the latter description, randomness is not the deciding 
factor in the resulting phase behavior. However, it is 
clear that our simulations are only compatible with the 
random-field hypothesis of de Gennes, since we see large 
differences between randomly and non-randomly placed 
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obstacles. 

Finally, having shown that domain formation in mem- 
branes with quenched obstacles belongs to the univer- 
sality class of the 2D-RFIM, further insights might be 
gained by comparing to other systems in that class. The 
obvious candidates are filled polymer blends confined to 
thin films [51], for which the connection to the random- 
field Ising model was already noted [52]. Experiments 
have revealed that quenched filler particles dramatically 
limit phase separation dynamics [53 . Complete phase 
separation is typically not observed, but rather a mosaic 
of microscopic domains [54]. Computer simulations and 
theory [52 , 55 , 56 indicate that the dynamics can become 



completely arrested, leading to many "pinned" domains 
[57l [58] . Some of these results are in remarkable agree- 
ment with the membrane simulations of Refs. [T8land[T9l 
and the present work. 
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